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摘 E 在 实际 工程 中 ， 由 有 限 元 模型 得 到 的 计算 值 与 通过 试验 获得 的 测量 值 之 间 往 往 存 在 偏 
差 ， 为 了 能 够 精确 预测 结构 的 动力 响应 ， 依 据 测量 信息 修正 现 有 的 动力 模型 是 非常 必 
要 的 。 本 文 研 究 结构 动力 模型 修正 中 的 一 类 对 称 和 矩阵 反问 题 (IP-MUP): 给 定 矩 阵 A = 
diag(w1,w2,-.-,w2) € R?*?, X = [z1, £2, ,z,] € R^*?, URIE Mo, Ko € SR^*^, K 
矩阵 M, K e SR"""^, {MXA = KX, X" MX —I,, HÆ M([Lr]) = Mo, K([1,r]) = 
Ko, K'PM([Lr), K([1,7]) SIIRRE M, K 的 前 > 阶 主子 矩阵 。 运 用 代数 特征 值 反问 题 的 
理论 和 方法 ， 文 中 给 出 了 问题 IP-MUP 有 解 的 充分 必要 条 件 ; 并 在 有 解 的 情况 下 ， 给 出 了 通 解 


的 显 式 表示 。 
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分 类 号 : AMS(2000) 65F18 中 图 分 类 号 : 0327; TB123 文献 标识 码 : A 
1 引言 
对 一 线性 无 阻尼 自由 度 系统 ， 其 动力 特性 可 由 如 下 常 系数 微分 方程 来 描述 


Mi(t) + Kaz(t) = 0, (1) 


其 中 z(t), 2(t) 413 n EAB ELE s EE LE. Ma, Ka 分 别 是 n x n BEAT Wr a St E RE 15 
刚度 矩阵 。 

在 实际 工程 中 ， 利 用 有 限 元 模型 (FEM) XX (1) 得 到 的 理论 分 析 结 果 与 实际 结构 的 实测 结果 
往往 存在 偏差 ， 这 种 偏差 主要 是 由 有 限 元 模型 误差 包括 边界 条 件 和 连接 条 件 的 简化 ， 几 何 模 
型 和 本 构 关系 的 不 准确 等 ) 以 及 实验 测量 误差 引起 的 ， 并 且 很 难 通过 改进 计算 方法 而 消除 凹 。 
近 二 十 年 来 ， 随 着 建 模 技术 和 动力 测试 技术 的 日 益 成 熟 ， 结 构 动 力 模型 修正 受到 了 科技 工作 者 
和 工程 师 们 的 广泛 关注 23。 修 正 的 结构 动力 模型 不 仅 能 够 更 加 精确 地 预测 结构 的 动力 响应 ， 
而 且 还 可 以 通过 结合 实测 结果 对 结构 进行 损伤 检测 和 剩余 寿命 评估 [7]。 

结构 动力 模型 修正 即 为 利用 结构 现场 实测 的 振动 信息 修正 结构 有 限 元 模型 ， 使 得 修正 后 
结构 分 析 的 模 态 参 数 与 试验 值 趋 于 一 致 。 假 定 由 模 态 参数 识别 技术 从 振动 测试 数据 中 得 到 了 
前 p 阶 精确 的 模 态 参数 ， 记 为 


A = diag{w?, w3, p) ER?*?, X = [z1, £2, +++ , £p] € RP, 
其 中 zi 分 别 是 第 ;i 阶 测量 频率 与 测量 振 型 ， 且 rank(X) = p。 现 在 要 用 模 态 参数 识别 结果 修 


正 初始 有 限 元 模型 ， 使 得 修正 后 的 质量 矩阵 M c SR” MARERE K c SR ”满足 下 列 模 
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态 正 交 条 件 
XTMX-—I, XTKX-A, (2) 


以 及 特征 方程 
MXA= KX. (3) 


鉴于 模型 误差 主要 来 自 结构 的 几何 形状 、 边 界 条件 和 受 力 状态 等 情况 复杂 部 位 ， 因 而 结构 
动力 学 模型 的 物理 参数 矩阵 中 常常 仅 部 分 元 素 存 在 明显 误差 ， 若 按 常 规 对 整个 模型 进行 修正 ， 
由 于 无 误差 的 元 素 也 参与 运算 将 会 导致 振 型 误差 的 引入 并 放大 ， 从 而 大 大 降低 修正 精度 。 为 了 
减少 这 种 影响 ， 并 降低 模型 修正 的 计算 工作 量 ， 只 需 对 存在 误差 的 元 素 进行 修正 。 借 助 于 误差 
定位 技术 8- 蔬 ， 可 确定 有 限 元 质量 矩阵 与 刚度 矩阵 中 存在 误差 的 元 素 。 通 过 重新 排列 有 限 元 模 
型 的 结 点 次 序 ， 可 使 质量 矩阵 与 刚度 矩阵 中 无 误差 的 元 素 集 中 在 和 矩阵 左上 角 的 一 个 区 域内 ， 
然后 利用 模 态 正 交 关 系 与 特征 方程 修正 误差 项 。 注 意 到 式 (2) 的 第 二 式 可 由 式 (2) 的 第 一 式 与 
式 (3) 推出 ， 由 此 ， 对 质量 矩阵 与 刚度 矩阵 的 同时 修正 即 为 如 下 的 矩阵 反问 题 。 

问题 IP-MUP 给 定 测 量 的 频率 组 成 的 矩阵 A = diag{w?, w ,ww2} € RP, BAW 
S WI d Æ ARE X = [21,22,77- , Tp] € R?*?, 以 及 矩阵 Mo, Ko € SR". KEREM, Ke 
SR"*"， 使 得 


MXA= KX, X'MX-1l, st. M([Lr]) = Mo, K([Lr]) = Ko, 


其 中 M([Lr]), K(Lr]) 分 别 表 示 和 矩阵 M, K 的 前 7? 阶 主子 和 矩阵。 

文 [12,13] 考虑 了 一 类 子 矩 阵 约束 下 的 矩阵 反问 题 ， 可 其 均 不 能 保证 修正 后 的 矩阵 满足 特征 
方程 。 文 [14] 考虑 了 一 类 振 型 矩阵 与 质量 矩阵 的 修正 问题 ， 但 其 并 没有 考虑 刚度 矩阵 的 修正 。 
本 文 运用 代数 特征 值 反 问题 的 理论 和 方法 ， 给 出 了 问题 IP-MUP 有 解 的 充分 必要 条 件 ， 并 在 有 
解 时 ， 给 出 了 通 解 的 显 式 表示 。 

ACHI Rnxm 表示 所 有 n x m 阶 实 和 矩阵 的 集合 ，SR"*x"* 表示 所 有 n x 即 阶 实 对 称 和 矩阵 的 集 
R, ORO RRN xn 阶 正 交 算 阵 的 集合 ，1 dos m Bris (a kE kE, AT 表示 和 矩阵 A 的 Moore- 
Penrose J” X if, 


2 问题 IP-MUP 的 解 
作 矩 阵 和 的 QR- 分 解 ， 得 


R 
x-o|1]. (4) 


Arb En x n MEXER, R e R” 是 非 奇 异 上 三 角 和 矩阵 。 由 XTMX = I,, & 


| RT, s]ewo| t| - (5) 
0 
id 
Mj M» p 
TMQ = , 6 
vues [8] n r 
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则 式 (5) 等 价 于 
M, M R 
lane) us «ME (T) 
比较 上 式 的 左右 两 边 ， 得 
Mi= R TR. (8) 
再 由 KXX = MXA， 得 
QTKQ | l = QTMQ | | A, (9) 
0 0 
id 
QTKQ = K, K> p 
T Ks n—p » (10) 


K, Ko R Mi M RA m 
KT K 0| | MZ Ms gs" 
比较 上 式 左右 两 边 ， 得 


Kı = M, RAR! = RTAR, Ki = MID, (12) 


其 中 万 = RAR-1。 从 而 可 解 得 


| RTR! M | 
M=Q QT, (13) 
MI M 
| R-TAR-! DTM: | 
K=Q QT, (14) 
MID Ks 
其 中 Mz € R?XQG-P 为 任意 矩阵 ，JM3, Ka € R(n-p)x(n—-p) 为 任意 对 称 和 矩阵 。 
由 和 M([1,7]) = Mo, K([1,r]) = Ko, 得 
E? 0 | M | z, o] =M, (15) 
| z, o| x | z, o] = xo. (16) 
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其 中 Qi e Rrxz，Q@a € Rrx(n-p), 将 式 (13), (14) 代入 式 (15), (16) 得 


R-TR- M QT 
p | u7 A jJ EON 
R-TAR-! DTM: T 
eeu me e 
即 
Q2M2 Qi + Q1 M2Q3 + Q2MsQ2 = Mo - QR T R 1QT Š H, (19) 
Q2 M2 DQI + Qi DT M;Q? +Q2K3Q3 = Ko —QiR TAR QI EG. — (20) 
VEERE Q> 的 奇异 值 分 解 


azul i (21) 


其 中 U = [Ui, Uj] e OR™, P = [P,, PB] e OR(-9*6-P, y; = diag(oi,-.-,0,), 8 = 
rank(Qo), o; > 0 (i = 1,--.,s), FE 


= Maj. Ma; 8 
P MsP -— , (22) 
MZ, M33 n—-p—s 
s n—p-s 
Ks! K32 s 
PTKP = (23) 
Ki, Kas n—p-s 


TEX (21) 与 式 (22) 代入 式 (19)， 得 


ZPLMIQTU, XPTMjQTUs UIQIM2P> 0 EMa£à 0 
十 * 
0 0 UIQuMSAX 0 


UTHU, UQHUs 
š : (24) 
UTHU, UTHU 
上 式 等 价 于 
Mai = X! (UT HU, — EP] MZ QIU: — UI Q1 MPE) D}, (25) 
UIQ MPX = UT HU,, (26) 


UZ HU = 0. (27) 
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同 理 ， 将 式 (21) 与 式 (23) 代入 式 (20)， 可 解 得 


Kasl = X (UIGU, - EPT MI DQTU, - UTQi DT Mj P.X)Y!, (28) 
UZ Q1DT M;P,X = UI GU,, (29) 
UT GU. = 0. (30) 


sÑ (26) 与 式 (29) 可 改写 为 


UTQ.,DT UTGU 
2 Qı MP = | ege (31) 
UF Qi UZ HU, 
记 
UTQIDT UTGU. 
geen sessio posee (32) 
UZ Q1 UZ HU 
并 设 和 矩阵 4 的 奇异 值 分 解 为 
arl š A (33) 
0 0 


EFT = [Tn, T] e ORC7722*0772, y = [Vi, Vj] e OR?**, A = diag(,---,0) t = 
rank(A), à; > 0 (i = 1,---,t). MEEDE (31) 有 解 Mo 的 充分 必要 条 件 为 中 


AATB= B (34) 
有 解 时 ， 其 解 M2 可 以 表示 为 
Ma = A* BP] + JPI + VoL, (35) 


其 中 Je Rox(-7-2, L e R(9-0x(-»n 是 任意 矩阵 。 

综 上 ， 已 经 证 明了 如 下 结论 。 

定理 1 给 定 测量 的 频率 组 成 的 矩阵 A = diag(u2,«2,---,w2) e Rzxp， 相 应 的 实测 振 型 矩 
阵 X = [z1, £2, £p] € Rnxp， 以 及 矩阵 Mo, Ko € SR7x7。 则 问题 IP-MUP 有 解 的 充分 必 
要 条 件 是 


UTHU;-0, UIGU,-0, AA'B-B, (36) 
HFG, H, Us, A, B r3 (19), 3X (20), =Ñ (21) 和 式 (32) 所 示 。 此 时 
R-TR- M. 
M-Q | : | QT, (37) 
MI M 


R TAR-1 DT M; T 
MID Ks 
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其 中 万 = RAR-1, Q WÑ (4) 所 示 
Ms = P Ms! Maj PT, K= P Ka! K32 pT, 
Mi, M33 KL K33 
Mai = X (UT HU, - PPE MI QIU: - UT Qi Mo B.E), 
Ksı = X (UT GU; - EP] Mj DQTU: - UT QI DT Ma P,>)> "1, 
P, M> 分 别 如 式 (21) 与 式 (35) MR, Mao, Ka» € R^* 79-9 为 任意 矩阵 ，Ma3，K33 € 
R(n-»-3)x(n-»-5 为 任意 对 称 和 矩阵 。 
3 XER 
例 ”描述 一 个 振动 系统 的 质量 和 矩阵 中 无 误差 元 素 组 成 的 矩阵 为 
52 22 18 -13 
22 
Mo — 0.03 x , 
18 13 104 0 


-13 — 0 24 


刚度 矩阵 中 无 误差 元 素 组 成 的 矩阵 为 


2 30 3 
Bg. 

Ko = 600 x 
-2 -3 4 0 


3 3 0 12 

测 得 系统 的 前 三 阶 固 有 频率 为 wi = 1.8841, wa = 10.0575, ws = 19.8203， 相 应 的 振 型 向 量 为 
T 
| 0.5585 —0.0841 0.3094 —0.0800 0.0996  —0.0553 0.0084 | ; 


T 
| 0.4751 —0.2353 —0.1717 —0.1646  —0.3562 0.0404 —0.1788 | ; 


T 
| —0.4241 0.2838 0.2512 0.0852 —0.0508 —0.2105 -0.4113 | ; 
容易 验证 给 定 的 数据 满足 定理 1 的 条 件 ， 据 式 (37) 与 式 (38) 求 得 问题 IP-MUP 的 一 个 解 M, K 38 
Æ |M XA — KX|| = 9.9695e — 013。 因 此 ， 测 量 模 态 数据 已 融 于 修正 模型 MXA—KX. 
4 结束语 


本 文 运用 代数 特征 值 反问 题 的 理论 和 方法 ， 研 究 了 结构 动力 模型 修正 中 的 一 类 对 称 矩 阵 反 
问题 。 数 值 算 例 表 明 ， 本 文 提 出 的 方法 不 仅 可 使 修正 模型 具有 对 称 的 结构 ， 而 且 由 于 避免 了 无 
误差 的 元 素 参 与 运算 ， 从 而 极 大 地 提高 了 结构 动力 模型 的 修正 精度 。 
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An Inverse Problem for Symmetric Matrices in Structural Dynamic 
Model Updating 
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Abstract: There is a gap between predictions resulting from finite-element model and experimental 
results of an actual structure. Updating the existing dynamic models based on the model test data 
is very important in order to predict actual behaviors of the structure. In this paper, an inverse 
problem for symmetric matrices in structural dynamic model updating (IP-MUP) is considered: given 
A = diag(w1,u2,::: v2) € R?'?, X = [z1,22,-- ,z,] € R^*?, and matrices Ko € SR°*", Mo € 
SR”*", find real n x n symmetric matrices K, M such that XTMX = Ip KX = MXA, with 
K([L1,r]) = Ko, M([1,r]) = Mo, where M([1,r]) and K([1,r]) are r x r leading principal submatrices 
of M and K, respectively. By means of the theory and method of the algebraic inverse eigenvalue 
problems, a sufficient and necessary condition under which IP-MUP is solvable is obtained, and an 
explicit expression for the general solution to the problem is also provided when the solvability condition 
is satisfied. 

Keywords: finite-element model; model updating; submatrix; inverse problem 


